Cosmography with the Sunyaev-Zeldovich effect and X-ray data 
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Cosmography provides a direct method to map the expansion history of the Universe in a model- 
independent way. Recently, different kinds of observations have been used in cosmographic anal- 
yses, such as SNe la and gamma ray bursts measurements, weak and strong lensing, cosmic mi- 
crowave background anisotropics, etc. In this work we examine the prospects for constraining cos- 
mographic parameters from current and future measurements of galaxy clusters distances based on 
their Sunyaev-Zeldovich effect (SZE) and X-ray observations. By assuming the current observational 
error distribution, we perform Monte Carlo simulations based on a well-behaved parameterization 
for the deceleration parameter to generate samples with different characteristics and study the im- 
provement on the determination of the cosmographic parameters from upcoming data. The influence 
of galaxy clusters (GC) morphologies on the Ho — qo plane is also investigated. 
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I. INTRODUCTION 



Cosmic acceleration is one of the most important is- 
sues of modern theoretical physics (see, e.g., Sahni & 
Starobinski 2000; Peebles & Ratra 2003; Padmanabhan 
2003; Alcaniz (2006) and Li et al. (2011) for recent re- 
views). However, more than a decade after its discovery, 
the physical mechanism behind this phenomenon remains 
unknown. In the context of Einstein's general relativity 
(GR), this result implies either the existence of a new 
field, the so-called dark energy, or that the matter con- 
tent of the universe is subject to dissipative processes 
(Lima & Alcaniz, 1999; Chimento et al. 2003). 

In order to investigate this phenomenon in a model- 
independent way cosmographic approaches have been 
successfully applied (see, e.g., Turner & Riess 2002; 
Bamba et al. 2012). This approach is independent on 
the gravity theory and the matter-energy contents fill- 
ing the Universe and provides an interesting way to map 
and explore the expansion history of the universe. Be- 
sides SNe la observations, other cosmological observables 
have been used to employ a kinematic description of the 
Universe as, for instance, cluster strong lenses (D'Aloisio 
& Natarajan 2011), SNe la plus cosmic microwave back- 
ground (CMB) and Baryon Acoustic Oscillation (BAO) 
(Santos et al. 2011; Demianski et al. 2012), BAO plus 
observational Hubble data and gamma ray bursts (Xu & 
Wang 2011), weak gravitational lensing plus gamma ray 
bursts (Wang & Dai 2011). 

On the other hand, the promising technique of mea- 
suring galaxy clusters distances based on their Sunyaev- 
Zeldovich effect (SZE) and X-ray observations have not 
been fully explored in the kinematic context. Our goal 
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in this paper is to derive cosmographic bounds on the 
Hq — qo plane by using angular diameter distance (ADD) 
from galaxy clusters. To this end we use a parametric ap- 
proximation of the deceleration parameter along the cos- 
mic evolution, given by q{z) — qo+qiz/(l+z). In order to 
ensure a period of structure formation during the matter- 
dominated era, we also make use of the asymptotic value 
of q(z) at high redshift, which reduces the above parame- 
terization to an one-parameter function. Initially, we use 
25 ADD of galaxy clusters as compilled by De Filippis 
et al. (2005). Since these authors obtained the ADD by 
using two morphological description for galaxy clusters 
(elliptical and spherical j3 models) we also explore the 
influence of the morphology on the results. By assum- 
ing the observational error distribution of the De Filippis 
et al. sample, we also perform Monte Carlo simulations 
to generate samples with different sizes and study the 
expected improvement on the determination of the cos- 
mographic plane Hq — qo from upcoming data. In what 
follows, we outline the main assumptions for our analysis 
and discuss the main results. 



II. GALAXY CLUSTERS DISTANCES 

The so-called SZE is a small distortion of the Cos- 
mic Microwave Background (CMB) spectrum provoked 
by the inverse Compton scattering of the CMB photons 
passing through a population of hot electrons (Sunyaev 
& Zel'dovich 1972). The effect is proportional to the 
electron pressure integrated along the line of sight, i.e., 
to the first power of the plasma density. Another im- 
portant physical phenomenon occurring in the intraclus- 
ter medium is the X-ray emission that occurs primar- 
ily through thermal bremsstrahlung. The X-ray surface 
brightness Sx is proportional to the integral along the 
line of sight of the square of the electron density. Briefly, 
in the context of these two phenomena one may consider 
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the different electronic density dependencies and evaluate 
the angular diameter distance of the galaxy cluster such 
that (Silk & White 1978; Cavaliere & Fusco-Fermiano 
1979), 



V{z) oc 
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(1) 



where Sxo is the central X-ray surface brightness, T e o is 
the central temperature of the intracluster medium, A e no 
is the central X-ray cooling function of the intracluster 
medium. ATq is the central decrement temperature, and 
9 C refers to a characteristic scale of the cluster along the 
line of sight whose exact meaning depends on the assump- 
tions adopted to describe the galaxy cluster morphology 
(Carlstrom, Holder & Reese 2002). This technique of 
measuring distances is completely independent of other 
techniques and it can be used to measure distances at 
high redshifts directly. Recently this technique has been 
applied to a fairly large number of clusters (Reese et al. 
2002; Jones et al. 2005; De Filippis et al. 2005; Bona- 
mente et al. 2006) - sec also Mason et al. 2001; Cunha, 
Marassi & Lima 2007; Lima, Holanda & Cunha 2010; 
Holanda, Cunha & Lima 2012, Holanda, Cunha, Marassi 
& Lima 2012 for current cosmological constraints using 
ADD from galaxy clusters observations. 

It worth mentioning that, in order to estimate the an- 
gular distance of a galaxy cluster from its SZE/X-ray 
observations, one needs to add some complementary as- 
sumptions about its geometry. In the last decade, many 
studies about the intracluster gas and dark matter dis- 
tribution in galaxy clusters have been limited to the 
standard spherical geometry (Reiprich & Boringer 2002; 
Bonamente et al. 2006; Shang, Haiman & Verdi 2009). 
The importance of the intrinsic geometry of the cluster 
has been emphasized by many authors (Fox & Pen 2002; 
Jing & Suto 2002; Plionis, Basilakos & Ragone-Figueroa 
2006; Sereno et al. 2006) and the standard spherical 
geometry has been severely questioned, since Chandra 
and XMM-Newton observations have shown that clus- 
ters usually exhibit an elliptical surface brightness. In 
recent papers, the elliptical morphology was also shown 
to be in agreement with the validity of the so-called cos- 
mic distance duality relation Dl = (1 + z) 2 Da, where 
Dl is the luminosity distance (Holanda, Lima & Ribeiro 
2010, 2011, 2012; Nair, Jhingan, & Jain 2011 ). 



III. COSMOGRAPHY 

Let us now assume that the Universe is spatially flat 
as suggested by current WMAP measurements (Hinshaw 
et al. 2012). In this case, the ADD in the FRW metric 
is defined by (c = 1), 



D A = (1 
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FIG. 1: 25 angular diameter distances from galaxy clusters. 
The red filled circle and the blue filled square correspond to 
elliptical and spherical description for the same galaxy clus- 
ters sample. 



where H(z) — a /a is the Hubble parameter and the de- 
celeration parameter q(z) is defined by 



ad dH 1 (z) 
^ ~d? = dt 



(3) 



In order to proceed further, we adopt in our analysis 
the following parametrization for the deceleration param- 
eter: 



q(z) = q + q 1 z/(l + z), 



(4) 



where qo is the current value of the deceleration parame- 
ter and in the infinite past q{oo) = qo + qi- This paramet- 
ric form was inspired by one of the most popular parame- 
terizations of the dark energy equation of state (Chaval- 
lier & Polarski 2001; Linder 2003) and, although very 
simple, seems to be flexible enough to mimic the q(z) be- 
havior for a wide class of accelerating models. However, 
for most of the viable cosmological scenarios, we expect 
the Universe to be matter-dominated at early times (i.e., 
after the radiation dominance), which implies q = 1/2. 
Thus, in order to ensure this past dark matter-dominated 
epoch whose existence is fundamental for the structure 
formation process to take place, we assume the constraint 
q{z » 1) = 1/2 (Santos, Carvalho & Alcaniz 2011). In 
this case, qo + qi = 1/2 and Eq. (3) becomes 



(5) 



which has the advantage of reducing our analysis to a 
two-parameter fitting, go and Hq. In what follows, we 
perform our analysis using this latter expression. 
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FIG. 2: Confidence contours on the plane h — qo for three data samples of ADD measurements of galaxy clusters and different 
morphologies. Left) 25 observational data of De Filippis et al. (2005). Middle) 100 simulated data points and Right) 500 
simulated data points. 



IV. GALAXY CLUSTERS OBSERVATIONS 

A. Current samples 

In order to derive cosmographic bounds on the epoch 
of cosmic acceleration we use the angular diameter dis- 
tances of galaxy clusters from the De Filippis et al. 
(2005) compilation. This sample is composed of 25 data 
in redshift range 0.023 < z < 0.8. These authors re- 
analyzed archival X-ray data with XMM-Newton and 
Chandra satellites of two samples for which combined X- 
ray and SZE analysis have already been reported using 
an isothermal spherical /3-model. One of the samples, 
compiled previously by Reese et al. (2002), is a selec- 
tion of 18 galaxy clusters distributed over the redshift 
interval 0.14 < z < 0.8. The another one (Mason et al. 
2001) has 7 clusters from the X-ray limited flux sample 
of Ebeling et al. (1996). De Filippis et al. (2005) used an 
isothermal elliptical /?-model and an isothermal spherical 
/? model to obtain Da(z) measurements for these clus- 
ters samples. As discussed in this reference, the choice of 
circular rather than elliptical j3 model does not affect the 
resulting central surface brightness or Sunyaev-Zeldovich 
decrement and the slope (3 differs slightly between the 
two models. However, significantly different values for 
core radius are obtained. The two samples of De Filippis 
et al. (2005) are shown in Fig. 1. The blue squares and 
filled red circles with the associated error bars stand for 
the same galaxy clusters but described by elliptical and 
spherical models, respectively. 



B. Simulations 

We have run a series of Monte Carlo simulations to 
generate synthetic samples of Da(z) for both the spher- 
ical and elliptical model. As a fiducial model, we used 



the best-fit values of qo and Hq obtained from statis- 
tical analysis with the current observational data, i.e., 
qo = — 0.85±i;|| and H = 75±10 km.s^.Mpc- 1 (ellipti- 
cal morphology) and qo = —X.17~ti'to an d Hq = 67.0±13 
km.s~ 1 .Mpc~ 1 (spherical morphology). Using these val- 
ues, we generated a number of numerically simulated 
ADD measurements. In order to estimate the error bars 
of the distance values, as well as the dispersion of the 
errors, we made a detailed study of the current data. 
Because of the relatively small number of observed clus- 
ters, a bootstrap analysis was carried out to make the 
numerical simulations more reliable. We assumed that 
at a given redshift, the values of Da are normally dis- 
tributed around the fiducial value, the error bars being 
drawn from a normal distribution as well. Samples con- 
taining 100 and 500 clusters were then generated. The 
statistical study of these synthetic data was then carried 
out and is presented in next section. 



V. ANALYSIS AND RESULTS 

In what follows, we perform a % 2 fit over the (qo, H ) 
plane. In our analysis we use a maximum likelihood that 
can be determined by a x 2 statistics, 



x 2 Wp) = E 



(D A (Zi;p) - V Ao ,; 



(6) 



'Vao 



where T> ao.i is the observational ADD, er-p Ao , is the sta- 
tistical uncertainty in the individual distance, Da is ADD 
as given by Eq. (2) and the complete set of parameters 
is given by p = (q , H ). 

The most common statistical error contributions to the 
ADD of galaxy clusters are: SZE point sources ±8%, 
X-ray background ±2%, Galactic N H < ±1%, ±15% 
for cluster asphericity, ±8% kinetic SZ and for CMBR 
anisotropy < ±2%. On the other hand, the estimates for 
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FIG. 3: Left: Normalized likelihood for qo from the three samples discussed in the text marginalizing over the Hubble parameter 
and considering a elliptical cluster morphology. Right: The same as in the previous panel for spherical morphology. 
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FIG. 4: Left: 2a bounds on the evolution of the deceleration parameter from observed (dot-dashed) and simulated (solid) 
samples assuming a elliptical cluster morphology. Right: The same as in the previous panel for spherical morphology. 



systematic effects are as follow: SZ calibration ±8%, X- 
ray flux calibration ±5%, radio halos +3% and X-ray 
temperature calibration ±7.5%. As a matter of fact, 
one may show that typical statistical errors amount for 
nearly 20%, in agreement with other analyses (Mason 
et al. 2001; Reese et al. 2002, 2004, see also table 3 in 
Bonamente et al. 2006). 

TABLE I: lcr estimates on the cosmographic parameters 



Data points 


cq /qo 


O HQ /Ho 


elliptical morphology 


25 (observed) 


1.58 


0.14 


100 (simulated) 


0.38 


0.08 


500 (simulated) 


0.31 


0.04 


spherical morphology 


25 (observed) 


1.15 


0.19 


100 (simulated) 


0.60 


0.13 


500 (simulated) 


0.42 


0.05 



In Figs. 2a - 2c we show the constraints on the plane 



(go, Hq) by using exclusively the SZE and X-ray observa- 
tions. In Panel 2a, we show the 68.3%, 95.4% and 99.7% 
contours for the De Filippis et al. (2005) samples with 
the solid and dashed lines corresponding to elliptical and 
spherical models, respectively. For these data and irre- 
spective of the model adopted, we find very loose bounds 
on the cosmographic parameters qo and Hq with both 
accelerating (q < 0) and decelerating (q > 0) universes 
being compatible at high confidence level. 

Our results show that, even considering the current 
observational error distribution, an increase in the num- 
ber of data points decreases considerably the uncertainty 
on the parametric plane go — Ho (see Table I) . It is also 
worth mentioning that the constraints on qo are insensi- 
tive to the galaxy clusters morphology, with the analyses 
with 100 and 500 data points ruling out a decelerating 
universe at 3a level. 

To quantify the improvement on the determination of 
the cosmographic parameters, we define our figure-of- 
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merit as FoM oc \ ja qQ uu a and find: 

F0M500 = 3.91 x F0M100 = 23.13 x FoM 25 (Elliptical) , 

F0M500 = 7.56 x F0M100 = 25.92 x FoM 25 (Spherical) . 

In Figs. 3a and 3b we plot the qo likelihoods for the 
two morphological descriptions when the Hubble param- 
eter is marginalized over. Clearly, a larger sample (100 
and 500 data points) with the same characteristics of the 
currently observed one can rule out a current deceler- 
ated universe. For the sake of completeness, we also plot 
in Figs. 4a and 4b a 2a reconstruction of the q(z) func- 
tion using both the current observational data (25 points) 
and 500 simulated data points. Differently from the cur- 
rent observations, the 500 simulated data points shows a 
switch from a decelerated to an accelerating phase. 

VI. CONCLUSIONS 

In this work we have discussed cosmographic bounds 
from measurements of the angular diameter distance of 
galaxy clusters based on their Sunyaev-Zeldovich effect 
and X-ray observations. By using a parametric approx- 
imation of the deceleration parameter along the cosmic 
evolution given by Equations (2) and (5), we have also 



explored the influence of galaxy clusters morphology on 
the estimates of Ho and qo. 

Two different analyses have been performed. First, 
we have derived bounds on the plane Ho — qo from the 
current sample of 25 ADD measurements for which the 
Sunyaev-Zel'dovich effect and X-ray surface brightness 
maps were fitted by an isothermal spherical /? model and 
an isothermal elliptical /? model. Irrespective of the mor- 
phology adopted, we have found very loose bounds on 
the space Ho — qo, with both accelerating (qo < 0) and 
decelerating (qo > 0) universes being compatible at high 
confidence level. In addition, we have performed Monte 
Carlo simulations of the ADD of galaxy clusters based 
on the current observational error distribution (De Fil- 
ippis et al. 2005) to study the dependence of the cos- 
mographic bounds with the size of the sample. We have 
shown that, even keeping the current observational un- 
certainty, an increase in the number of data points in- 
creases considerably the figure-of merit for the cosmo- 
graphic plane H n — qo- These results, therefore, highlight 
the cosmological interest in ADD measurements of galaxy 
clusters from the Sunyaev-Zeldovich effect and X-ray ob- 
servations and show that upcoming measurements may 
become complementary or even competitive with other 
cosmological probes. 
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